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We use computer simulations and a simple free energy model to study the response of a bi¬ 
layer membrane to the application of a negative (compressive) mechanical tension. Such a tension 
destabilizes the long wavelength undulation modes of giant vesicles, but it can be sustained when 
small membranes and vesicles are considered. Our negative tension simulation results reveal two 
regimes - (i) a weak negative tension regime characterized by stretching-dominated elasticity, and 
(ii) a strong negative tension regime featuring bending-dominated elastic behavior. This resembles 
the findings of the classic Evans and Rawicz micropipette aspiration experiment in giant unilamellar 
vesicles (GUVs) [Phys, Rev. Lett. 64, 2094 (1990)]. However, while in GUVs the crossover between 
the two elasticity regimes occurs at a small positive surface tension, in smaller membranes it takes 
place at a moderate negative tension. Another interesting observation concerning the response of a 
small membrane to negative surface tension is related to the relationship between the mechanical 
and fluctuation tensions, which are equal to each other for non-negative values. When the tension 
decreases to negative values, the fluctuation tension 7 drops somewhat faster than r in the small 
negative tension regime, before it saturates (and becomes larger than r) for large negative tensions. 
The bending modulus exhibit an “opposite” trend. It remains almost unchanged in the stretching- 
dominated elastic regime, and decreases in the bending-dominated regime. Both the amplitudes of 
the thermal height undulations and the projected area variations diverge at the onset of mechanical 
instability. 

PACS numbers: 


I. INTRODUCTION 


Bilayer membranes are quasi-two-dimensional (2D) 
fluid sheets formed by spontaneous self-assembly of 
lipid molecules in water [1], Their elasticity is tradi¬ 
tionally studied in the framework of the Helfrich effec¬ 
tive surface Hamiltonian for 2D manifolds with local 
principle curvatures ci and C 2 [2J 
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where the integration is carried over the whole sur¬ 
face of the membrane. The Helfrich Hamiltonian in¬ 
volves four parameters: the spontaneous curvature cq, 
the surface tension 00 , the bending modulus ko, and 
the saddle-splay modulus Rq. For symmetric bilayer 
membranes, Co = 0. If, in addition, the discussion 
is limited to deformations that preserve the topology 
of the membrane, then (by virtue of the Gauss-Bonnet 
theorem) the total energy associated with the last term 
is constant, and one arrives to the more simple form 


U= dS 
J A 


CO + -KO (Cl + C 2 y 


— CqA + —KqJ 2 , 


( 2 ) 

where A is the total area of the membrane and J, de¬ 
fined by J 2 = J dS (ci + C2) 2 , is the integrated total 
curvature. 


Helfrich Hamiltonian provides successful framework 
for describing many features of bilayer membranes and 


vesicles, including their large-scale shapes and trans¬ 
formations between them Q , membrane-membrane in¬ 
teractions Q , and membrane-mediated forces between 
proteins (“inclusions”) (5|. In light of this success and 
the wide acceptance off the model within the soft mat¬ 
ter and biophysics communities, it is surprising how 
poorly understood and confusing remains the surface 
tension term a 0 in Eq. ©■ Below we briefly review 
some of the complications associated with the concept 
of membrane surface tension. 

1. Membranes are much more flexible to bending 
than stretching. Therefore, in theoretical studies it is 
sometimes assumed that the lipids area density is fixed 
@, and changes in the total area result in from changes 
in the number of lipids In this picture, the surface 
tension is essentially the chemical potential of the sur¬ 
face lipids Q. Lipids, however, are highly insoluble in 
water (their critical micelle concentration is typically 
in the < 10~ 6 M range [lj), which makes the above 
interpretation for the surface tension largely irrelevant 
as there is almost no exchange of lipids between the bi¬ 
layer membrane and the embedding aqueous medium. 
More commonly, the surface tension is considered as a 
measure for the elastic response of a membrane with a 
fixed number of lipids to area variations 0,113 ■ This 
view clearly differs from the meaning of the term “sur¬ 
face tension” for fluid/fluid interfaces, where it serves 
as a measure for the free energy penalty resulting from 
an exchange of molecules between the “bulk” phases 
and the “interface” Q 

2. Further complication arises from the fact that the 
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membrane is a (quasi) two dimension manifold embed¬ 
ded in a three dimensional space. Therefore, it has two 
characteristic areas - (i) the total physical area A , and 
(ii) the area of its projection onto a planar reference 
frame, A p . Helfrich Hamiltonian involves an integral 
over the total area A ; but this area (whose determina¬ 
tion, taking into account the molecular structure of the 
membrane, is not without ambiguity) cannot be fixed 
due to the membrane thermal undulations. Thus, A 
is not a valid thermodynamic variable [l2j. It is, in 
fact, the projected area A p that emerges as the com¬ 
putationally relevant quantity B] . The force per unit 
length that needs to imposed on the frame in order to 
fix its area to A p is known as the mechanical (frame) 
tension, and will be denoted henceforth by t. Theo¬ 
retically, the mechanical tension can be identified with 
the derivative of the elastic free energy, F, with respect 
to A p : t = dF/dA p (where the differentiation is car¬ 
ried while holding the volume of the system, as well 
as the temperature and number of lipids, fixed). One 
can also fix r and let A p fluctuate. In this case, the 
relevant free energy is G = F — tA p , and the mean 
projected area satisfies (A p ) = —dG/dr. 

3. Another confusing aspect associated with the con¬ 
cept of surface tension is the distinction between Hel¬ 
frich Hamiltonian and the membrane free energy F. 
The latter is often referred to as the Helfrich free en¬ 
ergy , and the confusion arises because it is assumed to 
have the same form as the Helfrich Hamiltonian 

F = a A + t^kJ 2 . (3) 

In Eq. ©, A and J denote, respectively, the area 
and total curvature of the mean profile of the mem¬ 
brane. Eq. © features two renormalized coefficients, 
the tension a and bending modulus k, that are not 
equal to their intrinsic counterparts from Eq. © j er Q 
and kq. The statistical mechanics of thermal fluctua¬ 
tions around the mean profile are accounted for in the 
values of the renormalized elastic coefficients. For a 
membrane with a mean flat profile (i.e., not subjected 
to normal bending forces) A = A p and, J = 0. Then, 
Eq. © takes the form F = aA p , and r = 8F/dA p = a. 

4. Another quantity that can be identified as 
the membrane surface tension is the, so called, q 2 - 
coefficient 7 , also known as the fluctuation tension. 
The fluctuation tension is measured from the Fourier 
spectrum of the membrane height function with respect 
to the plane of the frame. For a membrane with a mean 
flat profile, the thermal average of the amplitude of a 
Fourier mode with wavevector q satisfies: (hf) = 0, 
and 

(\h4 2 \ = _ k eTA p _ 
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where fee is Boltzmann constant, T is the temperature, 
and l is a microscopic cutoff length. Some controversy 


surrounded Eq. © concerning the question whether 7 , 
the ^-coefficient in the denominator on the right hand 
side, is equal to a 0 or cr. This issue has been recently 
settled, and it now understood that the correct coeffi¬ 
cient is the renormalized surface tension a BIH , or, 
according to some theoretical studies, a slightly modi¬ 
fied version of this surface tension, ( A p /A)a B|. Sim¬ 
ilarly, the q A coefficient in Eq. © is equal to the renor¬ 
malized bending modulus k. Combining this with the 
discussion about Eq. ©. we conclude that both the 
measurable fluctuation and mechanical tensions coin¬ 
cide with each other (and with the renormalized ten¬ 
sion a): 7 = t = a. 

II. NEGATIVE SURFACE TENSION 

When the surface tension r = 7 = a vanishes, the 
membrane is “free to choose” the equilibrium projected 
area A p that minimizes the free energy G = F. The 
question we now wish to address concerns with the elas¬ 
tic response of the membrane to a further decrease in 
the frame area, which involves the application of a neg¬ 
ative surface tension. Based on Eq. ©, one may argue 
that for 7 < 0 , the membrane always becomes mechan¬ 
ically unstable because the amplitude of any mode with 
q < 1 J (— 7 /k) diverges. But such modes exist only in 
sufficiently large membranes, hence, small membranes 
can always sustain some negative surface tension. For 
instance, consider a square membrane of linear size L 
with k = 25 ksT ~ 10~ 19 J. From Eq. ©, one finds 
that such a membrane can withstand negative surface 
tension of size 7 = — 5 x 10 _ 3 N/m (which is comparable 
in magnitude to the typical positive rupture tension), 
provided that L < ( 27 t)a/— fc /7 ~ 30nm. This is the 
characteristic size of real small liposomes, and of bi¬ 
layers in highly coarse-grained simulations. Thus, the 
above estimation highlights the fact that the question 
of elastic response to negative surface tension is both 
interesting and relevant to current experimental and 
computational studies. 

The derivation of Eq. ©, and the proof that the 
fluctuation and mechanical tensions coincide with each 
other, involves several assumptions that do not neces¬ 
sarily remain valid when a becomes negative. Specif¬ 
ically, it is based on the investigation of the linear 
response of a mechanically stable flat membrane to 
small normal forces, and is restricted to configura¬ 
tions with smooth (twice differentiable) height func¬ 
tions h(f). But when cr < 0, the membrane can relieve 
the free energy cost of compression by buckling. We 
note that in the absence of normal forces (which is 
the case under consideration here), the system is not 
expected to undergo spontaneous symmetry breaking 
similar to that occurring e.g., in Ising spin systems be¬ 
low the critical point. The reason is that the membrane 
height profile h(r) is a continuous field and, therefore, 
the transition between different buckled configurations 
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(e.g., from buckled “upward” to “downward”) does not 
require the crossing of a free energy barrier. Thus, the 
system remains ergodic for negative tension, and due 
to the symmetry of the bilayer (h$) = 0 for all the 
Fourier modes. The questions that remain are: 

1. Does Eq. (U) still hold for cr < 0? We certainly do 
not expect it to remain valid for strongly compressed 
membranes since it is derived from a quadratic (in 
the height function) approximation to the full Helfrich 
Hamiltonian ©• However, considering the fact that 
it holds for a = 0 , there is no apparent reason why it 
should not hold for small negative a. 

2 . Are the mechanical and fluctuation tensions still 
equal to each other? (Obviously, this question is rel¬ 
evant only if the answer to question no. 1 is “yes”.) 
As noted above, the proof of this equality depends on 
the surface tension being positive. Now that it is neg¬ 
ative, the membrane prefers more buckled configura¬ 
tions with larger mean squared amplitudes. Does this 
imply that the fluctuation tension 7 drops faster (i.e., 
becomes more negative) than the mechanical tension 

T? 

3. What happens to the bending modulus k under 
compression? Recall that the coefficient appearing in 
Eq. (2J is the renormalized bending modulus which, 
just like the tension 7 , may vary with the frame area 
A p . For positive tensions, the variations in k are usu¬ 
ally negligible, but this may not be the case for neg¬ 
ative tensions when the membrane becomes increas¬ 
ingly more buckled. Does the increase in the degree 
of buckling under larger compressive stresses involve a 
decrease in k! 

4. Does the membrane exhibit linear (Hookean) elas¬ 
tic response to negative mechanical tension? In re¬ 
sponse to a positive tension, the membrane becomes 
stretched, and the relationship between the change 
in the area (strain) and the stress is indeed linear. 
However, the lipids constitute a dense two-dimensional 
fluid and, therefore, the membrane can be barely com¬ 
pressed below its most favorable physical area A$. 
When, under the application of a negative tension, the 
physical area A reaches A 0 , the negative tension causes 
the membrane to buckle and more and more area is 
“stored” in the out-of-plane fluctuations. This could 
lead to a highly non-linear elastic response. 

In what follows we will use Monte Carlo (MC) sim¬ 
ulations to address the above questions. 

III. COMPUTER SIMULATION METHODS 

To allow for large scale membrane simulations (N = 
1000 lipids per monolayer) in a computationally feasi¬ 
ble manner, we employ the Cooke-Deserno model 0 , 
with model parameters used in our recent work on 
fluid charged membranes [ 3 - I n this highly coarse¬ 
grained model, the lipids are represented as trimers 
with one hydrophilic and two hydrophobic beads, and 


the embedding solvent is handled implicitly via ef¬ 
fective interactions between the hydrophobic beads. 
The membranes are simulated in a box of linear size 
L x = L y = L, with periodic boundary conditions in the 
x — y (frame) plane. Initially, we place a pre-assembled 
flat membrane at the center of the simulation box, and 
we then allow it to equilibrate for 1 x 10 5 MC time 
units. On average, each MC time unit consists of N 
translation (with additional small intramolecular dis¬ 
placements) and N rotation move attempts carried on 
randomly chosen lipids. The membrane is simulated 
at constant frame tension r, which is accomplished by 
incorporating several move attempts, per time unit, to 
change the frame area of the membrane [3|. Each 
time unit also includes several collective “mode exci¬ 
tation” moves [ 2 OI] that accelerate the slow dynamics 
of the long-wavelength Fourier modes. For each value 
of r, the initial thermalization period is followed by 
a period of 1.5 x 10 6 MC time units during which we 
sample quantities of interest at 200 MC units intervals. 

In what follows, we set the thermal energy fcpT to 
be the elementary energy unit, and the length parame¬ 
ter of the Deserno-Cooke repulsive potential b to be 
the length unit. We simulate the membrane under 
frame tensions satisfying —0.3 < r < 0.5 (in ksT/b 2 
units) which, as revealed by our computational obser¬ 
vations, is the stability range of the simulated mem¬ 
branes. Relation to physical units can be made by 
setting b = 0.65 nm, which gives the unit of the sur¬ 
face tension /cpT/6 2 ~ 10 mN/n. For r > 0.5, the 
membranes rupture, while for r < —0.3, they ex¬ 
hibit large normal undulations leading to collapse and 
dissociation of lipids. Within the stability range, we 
measure the mean and variance of the projected area 
distribution ((A p ) and (6A 2 ) = (A 2 ) — (A p ) 2 , re¬ 
spectively). We also calculate the Fourier transform 
of the height undulations, by dividing the membrane 
into 8x8 grid cells, and calculating the local mean 
height of the bilayer within each grid cell. The Fourier 
transform of h{f) in wavenumber space n = q(L/2n) 
(ft = (n x , n y ) ; n x , n y = —4, —3,..., 2, 3) is defined by 

^ = iE^ e_2 ™' r7L - ( 5 ) 

r 

Notice that the linear size of the frame L 1 appearing 
(twice) in the definition of hfi, is not constant, but 
rather fluctuates during the course of the (constant 
tension) simulations. At each measurement, we use the 
instantaneous value of L. Also notice that defined 
in Eq. J5]) is dimensionless, due to the L _1 prefactor 
that does not exist in the more commonly used h$ of 
Eq. (|3|. In terms of the variable h^, Eq. dj) takes the 
form 

<l^| 2 >= ( 7 ) [ 7 (A P )(2^n) 2 + K (27m) 4 ]' (6) 
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FIG. 1: The spectral intensity as a function of the 
wavenumber for membranes under frame tension of r = 
—0.24 (squares), 0 (circles), and 0.24 (triangles). The 
dotted-dashed curves represent the best fits of the results 
to Eq. © over the first four modes. The inset shows the 
mean projected area per lipid as a function of r. 



There are four different modes corresponding to each 
value of |n|, and this number is reduced to two if 
In^l = \n y \, or if one of the components of n is zero. 
In the following section, the results for (and 

other related quantities) represent averages over these 
distinct modes. In Eq. ©, / is the grid size, which im¬ 
plies that L/l = 8 , independently of the instantaneous 
value of L. 


Due to molecular-scale protrusion, the physical area 
of the membrane cannot be unambiguously deter¬ 
mined. Therefore, we use the following approximation 
for (A) 


(A) 


(A p ) 




which is the physical area “visible'’ up to the resolu¬ 
tion of the grid. One can also define the effective area- 
stretch modulus of the membrane, Ka , by assuming 
that the free energy cost due to small variations in 
the projected area from (A p ) can be approximated by 
the quadratic form [1/2)Ka[A p — (A p )] 2 /{A p ) Un¬ 

der this approximation, the coefficient Ka can be ex¬ 
tracted from the fluctuation statistics of A p using the 
equipartition theorem 


_ k B T {A p ) 

m) 


( 8 ) 



FIG. 2: The fluctuation tension 7 as a function of the frame 
mechanical tension r. The solid line represents the equality 
7 = 7 - which is expected to hold for positive tensions. The 
inset shows the spectral intensity of the longest Fourier 
mode (n 2 = 1 ), (|/ii| 2 ), as a function of r. 


IV. RESULTS 


In section nu we raised several questions concern¬ 
ing the elastic and fluctuation behavior of membranes 
under negative mechanical tension. Here, we present 
computer simulations results addressing those ques¬ 
tion. 

We begin with the question of the validity of Eq. © 
for negative frame tensions. Figure [1] displays our re¬ 
sults for the fluctuation spectral intensity, \hft\ 2 , as 
a function of n 2 for membranes under three differ¬ 
ent mechanical tensions r = —0.24, 0, 0.24. The 
fits of the computational results to Eq. © are dis¬ 
played with dotted-dashed lines. Our analysis reveals 
that within most of the range of mechanical stability 
—0.3 < r < 0.5, the quality of each fit is very good. 
This demonstrates that Eq. © adequately describes 
the fluctuation behavior of bilayer membranes under 
both positive and negative tensions. 

From the fitting curves, one can extract the values 
of k and j{A p ) as a function of r, and by indepen¬ 
dently measuring the mean projected area, (A p ), one 
can obtain the value of the fluctuation tension 7 . Our 
attempts to use k as a single fitting parameters by forc¬ 
ing 7 = t resulted in poor fitting for negative tensions. 
The mean projected area as a function of r is plotted in 
the inset of Fig. [I] The observed increase in (A p ) with 
r is anticipated and will be discussed in detail later. 
Figure[2]depicts the fluctuation tension 7 as a function 
of r. The values reported in Fig. [5] are based on fitting 
analysis over the four longest fluctuation modes (small¬ 
est wavenumbers), and the error bars represent the in¬ 
tervals over which the fitting parameters, 7 and k, can 
be (mutually) varied while still producing reasonable 
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FIG. 3: The bending rigidity k as a function of the frame 
tension r. The dotted-dashed line is a guide to the eye. 


fits up to the accuracy of the computational results. 
For non-negative tensions, our results agree very well 
with the relationship r = 7 . As noted above in section 
ED there is no reason for this equality to remain valid 
for negative tensions, and our analysis summarized in 
Fig. [2] reveals that, indeed, 7/7 when the tensions 
are negative. Our results demonstrate that 7 < t and, 
as also argued above, it is likely that the more rapid 
decrease in 7 compared to r is related to the tendency 
of the membrane to form buckled configurations under 
negative tensions. The equality between 7 and r is re¬ 
gained for r ~ —0.15, and then the trend changes, and 
7 becomes larger than r. 

A closer inspection of the behavior of 7 depicted 
in Fig. [2] reveals that the 7 vs. r curve may be di¬ 
vided into three regimes: (i) a linear 7 = r regime 
for r > 0 , (ii) a non-linear regime where 7 < r < 0 
for mildly negative frame tensions, and (iii) a plateau 
regime (7 ~ const) for larger negative values of r. 
Saturation of the negative tension for strongly com¬ 
press membranes was previously observed ( 22 | , and we 
will return to the issue later in this section when we 
discuss our results for the physical area of the mem¬ 
branes (Fig. |4]). The fluctuation tension in Fig [2] is ex¬ 
tracted from Eq. ©, where it appears in the coefficient 
47 t 2 7 (A p ) of the n 2 term in the denominator. Naively, 
one may expect the saturation of the fluctuation ten¬ 
sion 7 to result in the leveling-off of the fluctuation 
spectral intensity |/i^| 2 . However, our computational 
results indicate that the amplitudes of the normal un¬ 
dulations continue to grow for decreasing values of r, 
as shown in the inset of Fig. [2] This apparent dis¬ 
crepancy can be only partially resolved by the trend 
in (A p ) whose value is reduced by about 10% in the 
plateau regime of 7 . The main factor explaining the 
increase in the undulation amplitude in the constant 
7 regime is the decrease in the effective bending mod¬ 



FIG. 4: Measured area as a function of the frame tension r. 
Solids circles and squares denote, respectively, the results 
for the frame and total areas. The former was measured 
directly from the simulation, while the latter was derived 
from the computed data for the spectral intensity, by using 
Eq.Q The open circle marks the optimal area of a flat ten¬ 
sionless membrane, Ao. The dotted-dashed line is a linear 
fit for the results for (A p ), while the horizontal dotted line 
marks Ao. All areas plotted in the figure are normalized 
per lipid. 

ulus k, the value of which is plotted in Fig. [3] We 
remind here that n is not a material but a thermo¬ 
dynamic quantity. For a tensionless membrane, the 
thermal undulations reduce (renormalize) the bending 
rigidity by Afc = —( 3 / 47 r)fcsTln(L/l), which is a small 
correction [23[. For r < 0, the amplitude of the fluc¬ 
tuations increase and, therefore, this correction term 
should become larger (in absolute value), which ex¬ 
plains the drop in the value of k seen in Fig. [3] In 
other words, just like the rapid decrease in 7 , reported 
above in Fig. [2] for membranes under negative tension, 
the reduction in n is also related to the increasing ther¬ 
mal roughness of the membrane, and the tendency of 
the membrane to form more buckled configurations. 

The results of Figs. [2] and [3] point to an interest¬ 
ing difference between the elastic coefficients 7 and n. 
The former decreases faster than r for small negative 
tensions and levels-off at large negative tensions. The 
latter exhibits “opposite” behavior and remains fairly 
constant in the small negative tension regime, and then 
decreases for strongly compressed membranes. The 
crossover between the regimes occurs at r ~ —0.15. 
Some light may be shed on these observations by the 
results of Fig. [I] depicting the mean projected and to¬ 
tal areas as a function of r. The results for the mean 
projected area, ( A p ), were measured directly from the 
simulations, while the data for the mean total area, 
(A), was calculated using Eq. J7}. For r > 0, we ob¬ 
serve a nearly linear dependence of both ( A p ) (see also 
the dotted-dashed line) and (A) on r. This behavior 
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FIG. 5: The stretch modulus I\a, measured from Eq. J5}. 
as a function of the frame tension r.The dotted-dashed line 
is a guide to the eye. 


agrees very well with the experimental results of Evans 
and Rawicz who also measured linear elastic response 
of giant unilamellar vesicles (GUVs) under positive me¬ 
chanical tension J 2 JI . We must point, however, to an 
important difference between the origins of linear elas¬ 
ticity in GUVs and small bilayer membranes. In both 
cases the linear elastic response is energetic in nature 
and dominated by the area elasticity of the membrane, 
while the entropy and bending energy of the height 
fluctuations play a secondary role in the response to 
stretching. In GUVs, this happens after the height 
fluctuations have been ironed by a very weak positive 
tension scaling inversely with A p . In small membranes, 
the height fluctuations are not damped and, in fact, 
the simulation results in Fig. [I] reveal that the excess 
area “stored” in the height fluctuations, (A) — (A p ), 
decreases only weakly with r. This implies that the en¬ 
tropy and bending energy of small membranes do not 
vanish (as in GUVs under tension), but simply exhibit 
relatively weak dependence on the frame tension (and, 
therefore, contribute weakly to the elastic response). 

In addition to the simulations of fluctuating mem¬ 
branes, we also simulated a tensionless (r = 0 ) flat 
membranes by running a MC code with moves allow¬ 
ing only local protrusions of lipids, but completely sup¬ 
pressing the longer scale bending modes (i.e., ensuring 
hft = 0 for all n). For a flat membrane A = A p . The 
measured area of the flat tensionless membrane, Aq, is 
denoted by the open circle and the horizontal dotted 
line in Fig. [I] This is the area that minimizes the elas¬ 
tic energy of the membrane. Fig.[I]provides an interest¬ 
ing interpretation for the weak and strong negative ten¬ 
sions regimes. The weak negative tension regime is es¬ 
sentially a continuation of the positive tension regime. 
The mean area of a tensionless fluctuating membrane 
is slightly larger than Aq, which implies that, in fact, 


FIG. 6: The bending, stretching, and total free energies (see 
definitions in the text) per grid cell, as a function of mean 
projected area (A p ) per lipid. The data for the total free 
energy has been shifted vertically by 0.3 for clarity. The 
vertical dotted-dashed line marks the measured projected 
area for r = 0. 


the membrane is stretched despite the negative me¬ 
chanical tension. Therefore, the area-dependent elas¬ 
tic energy continues to decrease with r into the weak 
negative tension regime. The strong negative tension 
regime begins when the total physical area reaches Aq. 
Since the membrane constitutes a dense fluid of lipids, 
it cannot be much further compressed, and in order to 
maintain the total area at Aq, more area must be ex¬ 
pelled into the height fluctuations. At this point, the 
clastic response becomes dominated by the height fluc¬ 
tuations bending elasticity and entropy, and we begin 
to observe a reduction in the effective bending modu¬ 
lus, k, instead of a reduction in the fluctuation tension 
7 . The saturation of the membrane physical area, and 
its correlation with that of the surface tension, was pre¬ 
viously reported [22J. Here, we demonstrate that this 
occurs when(A) reaches the value of Aq, which pro¬ 
vided an intuitive explanation to these observations. 
Notice that the rapid decrease in the projected area 
(A p ) with r in this regime is simultaneous to the in¬ 
crease in the projected area fluctuations. The resulting 
rapid decrease in the effective stretch modulus Ka, de¬ 
fined by Eq. (j%|). is plotted in Fig. [5] The vanishing of 
Ka signals the onset of mechanical instability [24}. 


V. DISCUSSION AND SUMMARY 

Based on our results, we identify two negative ten¬ 
sion regimes with the following features: (i) For weak 
negative tensions, the fluctuation tension 7 drops 
somewhat faster than the mechanical tension r, which 
is in contrast to the behavior observed for positive ten- 
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sions where 7 = r. In this regime, the membrane is still 
effectively stressed since the physical area (A) > A 0 . 
(ii) For strong negative tensions, the fluctuation ten¬ 
sion saturates, but then the effective bending rigidity 
begins to fall. Additionally, we also find that the to¬ 
tal membrane area in this regime reaches the optimal 
value of Ao, and does not continue to drop much [25j . 

In the spirit of Eq. © for the elastic free en¬ 
ergy of positively stressed membranes, we can ratio¬ 
nalize our observations for negatively stressed mem¬ 
branes by writing the free energy as the sum of two 
terms associated with stretching and bending. The 
former is given by the quadratic form A stre tch = 
( 1 / 2 ) Ka [(A) — Ao] 2 /Ao, while the latter may be eval¬ 
uated by Abend = (1/2) (\ h n\ 2 ^ 0]. These 

two contributions, and their sum are plotted in Fig. [G] 
as a function of (A p ). The free energies in Fig. [G] pro¬ 
vide insight into the computational results in this work. 
First, we notice that the total free energy attains a 
minimum at ( A p )/N ~ 1.32 (marked by the verti¬ 
cal dotted-dashed line), which is the mean projected 
area measured for r = 0. This must be the case since 
r = dF/dA p . Second, we notice that the bending free 
energy decreases with r, which is expected since, upon 
stretching, the thermal bending undulations tend to be 
suppressed. The stretching free energy increase with 
A p in the weak negative and positive tensions regimes, 
i.e., when (A) > A 0 . Under strong negative tensions, 
Astretch vanishes, which is associated with the observa¬ 


tion that the total physical area remains at the optimal 
value and does not change in this regime. 

The crossover from bending- to stretching- 
dominated membrane elasticity has also been observed 
in micropipette aspiration experiments of GUVs. 
There are, however, several key differences between 
the elastic behaviors of small and large membranes. 
Negative mechanical tension destabilize the large 
bending modes of giant membranes, while small 
membranes can withstand a (size-dependent) negative 
tension that may be comparable in magnitude to 
the positive rupture tension. In giant membranes, 
bending-dominated elasticity is limited to extremely 
small positive tensions that are typically two orders of 
magnitude smaller than the rupture tension. In small 
membranes, the crossover from bending-dominated 
to stretching-dominated elasticity is smoother and 
occurs at small negative tensions. In other words, the 
stretching-dominated elasticity regime extends into 
negative tensions, which stems from the fact that at 
zero tension, the membrane is still slightly stretched. 
Bending-dominated elasticity is observed at larger 
negative tensions. It is characterized by a decrease 
in the effective bending rigidity and stretch modulus 
of the membrane that ultimately lead to mechanical 
instability and membrane collapse. 
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